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Abstract. We present results of the search for coincident burst excitations over a 
24 hours long data set collected by AURIGA, EXPLORER, NAUTILUS and Virgo 
detectors during September 2005. The search of candidate triggers was performed 
independently on each of the data sets from single detectors. We looked for two- 
fold time coincidences between these candidates using an algorithm optimized for 
a given population of sources and we calculated the efficiency of detection through 
injections of tcmplated signal waveforms into the streams of data. To this purpose 
we have considered the case of signals shaped as damped sinusoids coming from the 
galactic center direction. In this framework our method targets an optimal balance 
between high efficiency and low false alarm rate, aiming at setting confidence intervals 
as stringent as possible in terms of the rate of the selected source models. 
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1. Introduction 

The network composed by the three gravitational wave (GW) resonant detectors 
AURIGA[H E!, EXPLORER and NAUTILUS [3, 4J and the interferometer Virgo[5] 
(hereafter called Virgo-bars network) is heterogeneous, as its single components differ 
for spectral sensitivity (see figure [T]) and antenna pattern. 

In the past years, various searches for GW signals have been independently 
performed by individual detectors or by networks of resonant bars (IGEC, [HI [7], ROG 
[8]) or interferometers [101 [11] . In the latter cases the networks were homogeneous: 
almost same antenna pattern (neglecting a small misalignment), similar (within a factor 
2) integrated sensitivity and roughly same observed frequency range (or detection 
bands). Therefore, a GW burst would produce approximately the same response in 
all the detectors of the network (notably, irrespective of direction and polarization of 
the source). In such cases the magnitude of observed signals can be compared directly. 

Previous burst searches among detectors with different spectral sensitivity and 
orientation were performed by the TAMA and LIGO Scientific Collaboration [12] among 
interferometers and by the AURIGA and LIGO Scientific Collaborations [T31IT4"] between 
interferometers and a resonant bar. 

The proposed network search strategy for the Virgo-bars data analysis takes as 
a starting point the IGEC coincidence search for burst GW events. This search was 
innovative with respect to previous searches as it preserved the detection efficiency by 
selecting the detectors which, time to time, had comparable directional sensitivity for 
sources located at a given sky position. In that case, however, there was no optimization 
on detection efficiency and the analysis relied on identical antenna patterns for the 
detectors. Instead, for the Virgo-bars network, it is necessary to further develop the 
idea included in the IGEC strategy. The detection efficiency will be determined by 
studying the software injections (Mock Data Challenge, in the following referred to as 
MDC) of a given collection of target waveforms. The approach attempted in this work is 
to use the efficiency computation both to tune the analysis parameters and to calibrate 
the final results, a step missing in ref. [7j. 

The results we present here are obtained in the simpler case of fixed time coincidence 
windows and two-fold coincidences among different detectors pairs. The coincident 
counts, divided by the detection efficiency and by the observation time, become then 
observed rates (or upper limits on rates) for that particular source population. The 
relevance of this study is methodological due to the short observation time, the 
uncertainty on the detectors calibrations and to some approximations in the production 
of the MDC. 

The paper is organized as follows: in section 2 we introduce the target GW signals 
and the source population we are dealing with. An overview of the exchanged data is 
presented in section 3. Section 4 presents the results obtained from software injections 
of GW signals into the data and the related estimates of the single detectors detection 
efficiencies and time errors. The coincident search strategy adopted in this work as 
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Figure 1. Typical spectral density of calibrated noise for the three resonant bar 
detectors during 2005 and for the Virgo interferometer in Sept 2005. 



well as the background estimation method and confidence of detection are described in 
section 5. The results and final remarks are presented in section 6 and 7 respectively. 
Finally, we report in Appendix A a summary of the pipeline main steps and in Appendix 
B a complete calculation of the energy budget associated to the injected signals. 

2. Target signals 

The class of transient GW signals is extremely large; moreover, such signals may be 
generated by a large variety of astrophysical sources. In this scenario, we have chosen to 
constrain the source population to the ensamble of waveforms that can be analytically 
described as a damped sinusoid (DS) with central frequency ranging within the bars 
bandwidths (850-950 Hz) and characterized by decaying times spanning at most a few 
tens of milliseconds. This choice is due to the different spectral densities of the various 
detectors in the collaboration (see figure [I]), so that the interesting signals for our specific 
network are the ones whose power is concentrated in the bars most sensitive frequency 
range. 

A typical damped sinusoid (DS) waveform is described by the following template: 

u(t) oc e" t/r cos(2tt f t + ip ) (1) 

where fo is the central frequency and r the damping time. These signals can be 
produced for instance by a ring-down phase following the merger of two black holes [T3] . 
Other sources whose emission can be modelled by are f-modes from neutron stars. 
The f-modes could produce a wave with variable frequency and damping time, which 
may sweep inside the observed frequency band [TO]. To make a realistic detection, the 



energy release should be about 10~ 3 — 10 -4 M Q for a galactic event (see Appendix B.l ). 

The astrophysical model for our source population considers elliptically polarized 
signals (as sources angular momenta should have random directions with respect to the 
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line of sight to the earth) incoming from the Galactic Center. 

3. Overview of exchanged data 

The exchanged data consists of event lists corresponding to 24 hours of data taking, 
starting from GPS time 810774700, or UTC time 14 Sep 2005 23:11:27. This choice 
corresponds to the longest scratches of continuous acquisition during the so called "C7" 
run of Virgo, when AURIGA, EXPLORER and NAUTILUS where in stable operation. 

Each group exchanged the triggers found above a chosen threshold by their 
respective burst event search algorithms. No further tuning of parameters and 
amplitudes has been done at this stage: a cut based on the magnitude of the events 
can be optimally set up afterwards based on the relative performances of the detectors 
at any given time. This selection reduces the number of background events without 
severely affecting the efficiency for a specific injection class. 

Before exchanging all the data, the time information has been offset by a secret 
time shift within each group. This was done in order to prevent any bias which might 
arise by looking at the zero-delay coincidence counts in the tuning phase of the analysis. 

It has to be noticed that the amplitudes may suffer from a systematic error due to 
the calibration uncertainty of each detector. This error is declared to be at most ~ 30% 
for Virgo, ~ 20% for EXPLORER and NAUTILUS, and ~ 10% for AURIGA. 

3.1. Event Trigger Generators 

From the Virgo side, Power Filter[T7] was the chosen Event Trigger Generator. Power 
Filter searches on whitened data for a power excess using different time analysis windows 
and different frequency bands and it uses as an indicator of the signal magnitude the 
(logarithmic) Signal to Noise Ratio (SNR). Events were exchanged at (logarithmic) 
SNR > 3.4. 

AURIGA group has successfully tested on its data an implementation of WaveBurst 
Event Trigger Generator, which is an excess power algorithm based on the wavelet 
transform developed by the LIGO Scientific Collaboration |18j : the exchange threshold 
was set at amplitude SNR > 4.5. For the NAUTILUS and EXPLORER detectors the 
Event Trigger Generator is an adaptive linear filter matched to the impulse response 
[19l|20]. The amplitude calibrated for the impulse response and the SNR were exchanged 
for each event and the exchange threshold was fixed at SNR > 3.8. 

3.2. Data quality: the correlograms 

The histogram of time difference among the outputs of the local Event Trigger 
Generators (or correlogram) with a given bin size is equivalent to counting coincidences 
with fixed window as a function of the time delay. The absence of cross-correlation is a 
useful hint that the coincidences at some delay may be considered as representative of 
the accidental coincidence probability at zero delay. 
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Figure 2. Spectral densities of correlograms and cross-correlograms for AURIGA, 
EXPLORER, NAUTILUS and Virgo. The spectra are flat within errors for all cases 
but for Virgo, where a 0.6Hz noise line is dominant (this is visible as periodic ripples 
in the Virgo auto-correlogram 



While single detectors, like Virgo, show some level of auto-correlation in the 
correlogram of the background events, on the other hand the cross-correlogram is flat 
as expected for random Poisson point processes. 

We found no evidence of modulation in the cross- correlation histogram up to ±400 s. 
Looking at the histograms in the Fourier domain (see figure [2j, the spectral density is 
flat for all detectors and detector couples, with the exception of Virgo, which presents 
a known strong peak at 0.6 Hz; this is the fundamental pendulum mode of the Virgo 
mirror suspension, whose excitation is only partially suppressed by the interferometer 
control loops. 



4. Results from software injections 

The injected GW signals consist of time series of sampled DS with /o within the sensitive 
frequency region for the resonant detectors, r spanning between 1 and 30 ms (see table Q 
and random elliptical polarization. 

The source location is chosen at the Galactic Center. For each detector, a 
specialized time series is produced including the time delays and the amplitude 
attenuation due to antenna pattern (see for example the amplitude modulation on 
figure [6]), using the SIESTA simulation software |21j . The simulated signals arrive at 
the Earth center approximately evenly spaced by 10 s (with a random jitter of ±0.5 s), 
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producing a set of 8640 injections over the 24 hrs observation time. 

In the following, when referring to the "amplitude" of the population, we mean the 
absolute h rss amplitude of the wave, i.e. the amplitude at the earth of the unprojected 
wave tensor (see Appendix B.l ). The generated waveform amplitudes for the coincidence 



analysis range between h rss = 1CT 20 Hz -1 / 2 and 10~ 18 Hz _1//2 in order to span the curve 
of efficiency vs h rss for all the detectors - see section 4.3. 

We acknowledge a coarse approximation in the algorithm that calculates the 
samples amplitudes for the MDC, resulting in an overall underestimation of the injected 
signal strength h rss with respect to the declared value. The effect of such approximation 
is fully negligible for the most energetic signals, while it may cause a spurious loss of 
efficiency at the lowest amplitudes and for the less sensitive detectors. Nevertheless, the 
reported results are conservative and the methodological relevance of this work is not 
affected. 



4-1. Time errors 

The timing error of all search algorithms is heavily dominated by systematic biases. 
This is typical of algorithms that are not matched to the particular signal one is looking 
for. For instance, for the Virgo detector, the Power Filter filter bandwidth is ~ 100 Hz 
in the narrowest channel, and the time of arrival is determined by the time when the 
signal reaches its maximum amplitude. Because of this the biases for 1 ms and 10 ms 
long DS are 0.8 ms and 3.6 ms respectively. 

For the AURIGA detector, the bandwidth is narrower, causing a larger distortion 
of the signal. Moreover the time associated to the event is computed as the baricenter of 
the signal profile above threshold. The amount of the bias depends strongly on the time 
duration and also on the central frequency of the signal. Altogether, for the durations 
< 30 ms the bias for DS ranges from ~ 3 ms to ~ 30 ms depending on the central 
frequency [22]. 

The linear filter matched to a delta signal is unbiased for wide band signals, such 
as 1 ms long DS, but it is more biased as the signal duration gets longer [23]. For the 
considered DS the worst time error of EXPLORER and NAUTILUS was of the order 
of 8 ms and 16 ms for damping time r = 10 ms and r = 30 ms respectively. 



Table 1. The parameter space of DS signals is described by their central frequency 
/o and their decay time r. The latter takes values spaced logarithmically by about a 
factor 3, while the frequency axis is sampled at a special subset of frequencies, which 
was chosen on the basis of the typical narrow-band power spectral densities of the bars. 



r (ms) 


fo (Hz) 


1 


914 


3 


882 946 


10 


866 898 930 


30 


866 874 906 930 938 
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Figure 3. Distribution of time differences of detected events in coincident injection 
of DS(914Hz;lms) (left), DS(930Hz;10ms) (center) and DS(930Hz;30ms) (right), at 
h rss = 1(T 18 Hz" 1 / 2 for the couple AURIGA- Virgo. 

In figure [3] we see for example that the maximum time difference between AURIGA 
and Virgo is < 30 ms. 

4-2. Distribution of amplitudes of accidental events and injected signals 

The single detector search algorithms provide different estimates of the magnitude of 
the signal. Although not directly comparable among different detectors, the event 
magnitudes provided by each detector show how much the population of injected 
waveforms stems from the noise distribution. An example can be seen in figure |4} 

4-3. Efficiency of detection 

The software injections have been used also to monitor the detection efficiency of the 
single detectors. The efficiency is computed for different waveform amplitudes and using 
a ±40 ms time window around the injection times: we look for coincidences between the 
injected signals and the events found by the Event Trigger Generators. The calculation 
of the efficiency is based on the nominal 24 hours allocated for this search. Consequently, 
dead times in the data due, for example, to epoch vetoes affect the average efficiency. 
This is the case for detector Virgo: its duty time is, in fact, far less than 100% in the 
24 hours considered (about 7 hours are vetoed out because of bad data quality), which 
finally reduces the attainable average efficiency (see figure [5J. The resonant detectors 
have instead a very stable duty cycle and they show improved performances with respect 
to Virgo when the signal is fully contained in their bandwidth (which requires special 
selection of central frequency and long signal duration). 

In conclusion, our definition of efficiency includes also dead times when data are 
missing or vetoed out from the nominal set, as we are interested in whether or not the 
network was able to recover the injected signals. We will comply with this comprehensive 
definition of efficiency throughout the paper. 



First joint Gravitational Waves search by the AURIGA-EXPLORER-NAUTILUS-Virgo collaborationlO 



-background 

-DS(914Hz.1ms) 

-DSl930Hz,10ms) 




12 16 20 24 
magnitude 



background 
DS(914Hz,1ms) 
-DS(930Hz,10ms) 




3 5 5 6 5 



9 5 11 12 5 14 15.5 17 18.5 20 
magnitude 



background 
DS(914Hz,1ms) 
-DS(930Hz,10ms) 





Figure 4. Distribution of detected "event magnitudes" for background events 
and injections of DS(914Hz;lms) and DS(930Hz;10ms) at h rss = lO" 19 ^" 1 / 2 for 
AURIGA (top-left), EXPLORER (top-right), NAUTILUS (bottom-left) and Virgo 
(bottom-right) . The "magnitude" can be the SNR ratio given by WaveBurst algorithm 
(AU), the one given by a linear matched filter (EX and NA) and the logarithmic SNR 
of Power Filter (Virgo). Given the very different meaning of these quantities, the plots 
obtained for different detectors cannot be compared directly. 




1.E-21 1.E-20 1.E-19 1.E-18 1.E-17 1.E-21 1.E-20 1.E-19 1.E-18 1.E-17 

waveform amplitude (hrss) waveform amplitude (hrss) 

Figure 5. Efficiency of detection for the four detectors when recovering injections of 
DS(914Hz;lms) (left) and DS(930Hz;10ms) (right) at different values of h rss . In the 
selected 24 hours, ~ 7 hours of Virgo data have been excluded by epoch vetoes based 
on data quality. That is why efficiency levels at about 70%. 
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5. Coincidence analysis 



A coincidence between two detectors is denned as the fulfilling of the relation 



(*) Ah) 



At 



< T 



(2) 



where t{!? is the estimated time of arrival of the jk signal in the detector labelled 



jk 



k and Athk is the time of flight between the sites of the two detectors. For the case of 
our network and for signals coming from the Galactic Center, the maximum time delay 
ranges between 0.3 ms for the couple Virgo-AURIGA and 2.4 ms for EXPLORER- 



NAUTILUS. From the results of the MDC injections (specifically, sec, 4.1), we set in the 
following T w = 40 ms. After the tuning of the analysis (which will be described in the 
following sections), we checked, for all the configurations shown on tabld^J that the loss 
in the overall efficiency with this coincidence time window is at most 1%. 



5.1. Thresholds optimization 

As we can see from figure |4j for a small increase in the magnitude threshold a large 
reduction of the background counts can be expected. The accidental coincidence rate 
between two detectors is proportional to the event rate of the two detectors, so we can 
act on one of the two thresholds, or on both. The trade off is the reduction of detection 
efficiency. Sometimes the detected magnitude of the injected events is large enough to 
increase the threshold up to exclude any background coincidence, while in other time 
periods lower thresholds are preferred to preserve the detection efficiency. 

In order to quantify these statements, we consider a gain function defined as the 
ratio between the average efficiency and the square root of the background counts. The 
rationale for this choice derives from the procedure to set confidence intervals on the 
number of true coincidences. In fact, the background of accidental coincidences can 
always be subtracted from the found (total) number of coincidences. The residual of 
the subtraction is -loosely speaking- the number of truly correlated events between the 
detectors. In this sense, our gain function is the ratio between the average efficiency 
and the fluctuations of the background. As our source population is set in the Galactic 
Center, we apply a time varying threshold (calculated every 30 minutes by maximizing 
the gain function, as above stated). This implies that for each time bin the threshold 
is set at the level corresponding to the maximum of the gain function for the pair of 
detectors. The overall result is that we apply over the entire data set and for each 
detector of the couple a non-constant cut on the event magnitude, using a threshold set 
every 30 minutes (the analysis pipeline is discussed in more detail in Appendix A). 

Figure [6] shows an example of such adaptive thresholds. As Virgo, in this example, 
can see the injected events much better than AURIGA, the algorithm starts with 
raising Virgo threshold. This cleans up most of the coincident events. The efficiency of 
AURIGA is therefore preserved, as its threshold is left almost untouched. 
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5.2. Background and efficiency of the network 

The efficiency of detection is empirically defined by the sets of data containing MDC 
injections. The ratio between the coincident events due to injections found in the 
detector couple under study and the known number of injected events gives the empirical 
estimate of the efficiency [|J 

For the background estimation, we first take care of the (possible) true correlated 
events present in the data by shifting the times of the event lists of different detectors 
before looking for coincidences in time. By repeating this operation a number of times, 
we get renewed instances of the counting experiment: we will refer to this procedure as 
time delay analysis. Altogether, the coincidences from hundreds of shifted configurations 
provide a rich population from which we can determine the main parameters of the 
background distribution. If the time slide measurements are independent from each 
other, the number of accidental coincidences in each time shift should be Poisson 
distributed. We tested this hypothesis by means of a \ 2 test just on those searches 
which have an high expected number of accidental coincidences, Nj,, so to ensure a 
sufficiently large data sample. The corresponding p-values were not inconsistent with 
the Poisson model for the expected number of accidentals. 

The optimization procedure described in the previous section determines the cuts 
on the data set based on a function of the estimated background and efficiency. Hence 
these two estimates will be biased, sometimes severely, and cannot be used for setting 
confidence intervals. For this reason, we preliminarily divide the original data (accidental 
coincidences and MDC injections) into two equal size subsets: one is used in the 
optimization phase. The threshold levels obtained at the end of the optimisation 
are subsequently applied to the second halves of the data, without further tuning. 
Background counts and efficiency are thus computed from this second subset, giving 
unbiased estimates. 

5.3. Setting confidence intervals 

We first set the confidence interval on the number of correlated events detected in 
coincidence, following a unified approach in the spirit of [21]. However, the procedure 
we adopt to build the confidence belt is different and its fundamentals has been discussed 
in [7] and [25] ■ We start by considering the likelihood of the number of coincidences 
at zero-delay N c as a function of the expected values of the accidental counts and of the 
correlated events, N b and N GW respectively: 




(3) 



The confidence intervals on Nqw are built by integrating the likelihood 




(4) 



| The contamination due to accidental coincidences was found to be negligible. 
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Figure 6. Example of threshold placing for detection optimization for the template 
DS(914Hz;lms) at h rss = 1CT 19 Hz" 1 / 2 for the couple AURIGA (above) and Virgo 
(below). The abscissa represents the period of the exchanged 24 hours of data taking, 
the ordinate is the event magnitude given respectively by the WaveBurst and the 
PowerFilter algorithms. The light markers indicate background coincident events, 
the darker ones are events generated by injections. The modulation of the event 
magnitude by the antenna pattern is clearly visible in the two plots. The events 
form coincident couples, i.e. to each event in AURIGA the corresponding event in 
Virgo is plotted. Excluding one event from one detector automatically excludes the 
paired one in the other. Thus, given the relatively good separation of MDC events 
from background events in Virgo, the threshold here is raised up to the limit where 
basically all background events are excluded, while the threshold in AURIGA is left 
almost at its initial value. 
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Figure 7. Measured global confidence of detection over all configurations vs nominal 
confidence of the null hypothesis test for each single search. The circles refer to the 
whole list of trials, the squares to the configurations which gave the best upper limits. 

to find the smallest interval [N in f, N sup ] corresponding to the chosen / value. The 
set of these intervals computed for the possible N c values makes the wanted confidence 
belt of Now vs N c for fixed level of background iVj,. It is well known that the / cannot 
be interpreted as a frequentist probability; therefore, the actual probability that the 
confidence intervals [N in f, N sup ] include the true New value, i.e. the coverage, must be 
empirically determined by a Monte Carlo [25]. The quoted coverage is the minimum 
coverage ensured by the belt [§J 

In the present work, we modify the above procedure by adding a more stringent 
test of the null hypothesis in order to increase the coverage when N GW = 0. In practice, 
we want a more stringent false alarm to issue a two-sided confidence interval with lower 
extreme greater than zero. This is done by performing a Poisson one-tail test on the 
found coincidences N c assuming the null hypothesys, New = 0. The significance of the 
test is set at a higher level than the coverage of the belt, as discussed in the next section. 
In case the test is passed, the lower extreme of the confidence interval N in f is extended 
to zero. 

§ Notice that at 95% coverage, the upper limit set by this belt is 3.6 counts when N c = 0, regardless 
of the value of the background, iV&. 
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5.4- Multiple trials and global confidence 

In this analysis we perform many searches for different signal waveforms and amplitudes 
by different detector pairs. The trial factor is therefore large and we are interested in 
controlling the global probability of false claim of the whole experiment. To estimate it, 
we considered the 400 different time-delayed configurations obtained by time shifting the 
original data set, assuming they are independent realizations of the experiment with no 
correlated events. The resulting coincidences has to be accidental and we can empirically 
estimate the distribution of these coincidence counts. We can then determine the global 
confidence of detection of the entire experiment as a function of the chosen significance 
of the null hypothesis test on each single trial. The results are shown as circle data 
points in figure [7} In particular, to have a global confidence of detection of at least 95%, 
we need to set the significance up to ~ 0.999 on the null hypothesis test for each single 
trial. 

There is a trade off between the global false claim probability and the detection 
efficiency of the experiment: to decrease the false claim probability it becomes more 
difficult to recognize a true signal. For instance, with ~ 0.1, as observed in most of 
the configurations, and the threshold of the 1-tail Poisson test set at 99.9%, we need at 
least 3 coincidences detected in a couple of detectors in order to reject the null hypothesis 
and make a claim for correlated events. 

In order to reduce this drawback, we attempted to limit as much as possible the 
trial factor of the experiment. In particular, for each signal waveform and amplitude we 
considered only the results produced by the pair of detectors which was performing best 
(square data points in figure [7]). To select the best performing pair, we compute fake 
upper limits per each pair under the assumption that the number of found coincidences 
is compatible to the expected background of the pair. Only the pair producing the more 
stringent fake upper limit is then searched for true coincidences. In this way, setting a 
significance of 0.999 on the null hypothesis test for each single trial, we achieve a global 
confidence of ~ 99%, i.e. a false claim probability of ~ 1%. 

6. Results 

As a last step, we searched for coincidences at the true time in all the selected 
configurations. The final outcome was consistent with no rejection of the null 
hyphothesis test at 99.9% confidence for each configuration, corresponding to a global 
false claim probability of 1%. The confidence interval are, hence, upper limits on the 
rate of incoming GWs. 

Detailed results are presented in table [2j where, for each injected waveform, the 
estimated efficiency, the average background, the zero-lag coincidence counts and the 
corresponding 95% upper limit are reported. Rescaling the listed efficiencies by their 
asymptotical values, we can easily infer that, for the 11 injected waveforms, the so- 
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Figure 8. Upper limit at 95% conservative coverage on the event rate as a function 
of the population amplitude. Upper limits referring to the same duration of the signal 
but different central frequencies are grouped with the same symbol. Notice that the 
systematic error of about 10-30% on the amplitude calibration has not been taken into 
account. 



called [jjj ranges from 5-10 20 Hz 1//2 (for the DS with central frequency /o=930 Hz 
and r=30 ms for the AURIGA-NAUTILUS pair) to 1 ■ 10~ 19 Hz" 1 / 2 (for the DS with 
central frequency /o=914 Hz and r=l ms for the AURIGA-Virgo pair). We included 
also a few configurations with low amplitude signals: in order to preserve the residual 
efficiencies, the related average backgrounds are quite high, leading to non null zero- 
delay coincidence counts. 

Figure [8] shows the upper limits as a function of signal amplitude. The asymptote 
for large amplitude signals is inversely proportional to the observation time and to the 
asymptotical efficiency and depends on the confidence belt: for this 24 hr search, it is 
~ 4.0 events/day, as our maximum efficiency is ~ 90% and the chosen confidence belt 
(see sec. 5.3[ ) with 95% C.L. sets a pretty conservative value (3.6 counts). 



|| The signal amplitude with 50% detection probability. 
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Table 2: Results of the two-fold coincidence searches, 



for the chosen couples of detectors (see sec 5.4) 
and for each set of the waveform parameters and 
amplitudes (Second, third and fourth column). The 
last four columns represent, respectively, the efficiency 
of detection measured by means of MDC injections, 
the average accidental coincidence counts (iV b ), the 
number of coincidences found at zero-delay (N c ) and the 
corresponding 95% upper limit. The errors associated 
with the efficiency and background estimates show the 
la statistic fluctuation (apices). 
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7. Final remarks 



We presented a methodological study for analyzing data collected by a network of non- 
homogeneous detectors. The search was aimed at detecting transient GW signals. We 
implemented a two-fold time coincidence search; however, this method could be applied 
as well to any detectors combination, e.g. three-fold, logical "OR" of two-fold, etc.. For 
each set of waveform parameters and amplitude, 6 different couples of detectors were 
available: we chose to perform our search on those couples which allow potentially the 
set up of the most stringent upper limit (see sect. 5.4). 

The key point of the method is the optimization process of the analysis thresholds 
for a given source population by means of Monte Carlo MDC injections. 

Although the proposed methodology is viable for any specific signal model, 
including the sky distribution of the sources, in the preset study we assumed DS signals 
incoming from the GC, limiting our observation range to our galaxy. 

Moreover, in order to estimate the detection efficiency of the network we applied a 
standard Monte Carlo procedure based on a large set of injected signals. Although this 
software technique is computational intensive, it permits to derive reliable values of the 
efficiency and unbiased physical interpretation of the results, i.e. the GW amplitudes 
and rates of the population under study. 

Finally we notice that in the procedure presented here, the statistical test relies just 
on the event magnitude. However, we stress that it is possible to extend the method 
either by including other statistical tests in the definition of the local Event Trigger 
Generators or by implementing a common maximum likelihood estimator (x 2 test). 
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Appendix A. Optimization pipeline 



We present here a schematic overview of all the steps: 
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• A specific search algorithm for each detector is run in order to produce lists of 
triggers. For each trigger the time and an estimate of the SNR ratio are exchanged. 
The exchanged lists comprise one derived by the analysis of the plain data from 
the detector, and others obtained by adding to the data different MDC channels 
before running the event search algorithm. 



Coincidence search 

• For each couple of detectors, the plain data sets are searched for coincidences after 
adding 800 offsets in steps of Is, with a time window of 40ms. This covers about 
7 minutes before and after the unshifted time, with a safety range of ±20s around 
the zero-delay tim^U 

• All the lists with MDC injections are searched for coincident events (obviously no 
time delay analysis is performed in these cases) and the detection efficiency is thus 
evaluated. 



Optimization procedure 

• The data are split into two equal sets, one used in the optimization phase, the other 
in the estimation phase, by alternating shift index. 

• The data are divided in 30 minutes long time bins . 

• We evolve a couple of staircase thresholds (jointly for the two detectors) by 
increasing the threshold level in one time bin and one detector at a time. At 
each step, in order to have a significant variation, the test threshold is increased by 
an amount which corresponds to a reduction in the background counts of the order 
of the standard deviation of the counts themselvea3 

• At the nth step of the algorithm, we compute the ratio N^l / \J N^ n \ where 

and are the total number of MDC coincidences and background coincidences 
from the set reserved for the optimization whose associated amplitudes are above 
the nth set of thresholds. Then the effect of increasing the threshold of one level at 
one bin is evaluated by computing the new ratio N^ l) / \J Nl n+l) . Every time the 
bins in both detectors are tried one by one in order to find for which a threshold 
change would score the higher benchmark. If this benchmark is better than the one 
obtained at the previous step, the level for that bin is changed, and this is taken as 
the starting point for the next loop. Instead, if all changes resulted in a decrease of 

% We recall that 10s is the maximum value for the blind shift applied previously to data exchange, see 
section [3j 

+ We found empirically that the square root of the counts divided by 6 is a good compromise. Yet, we 
impose that the background decrease of at least 4 counts. 



First joint Gravitational Waves search by the AURIGA-EXPLORER-NAUTILUS-Virgo collaboration21 



the benchmark, the loop is exited. The loop will otherwise continue up to reaching 
the higher level of thresholds. 

• When a time bin is found contributing to the efficiency by less than 0.1% the 
threshold in that bin is raised until no background event is surviving. Similarly 
to what was done, at the end of the previous loop we recompute the benchmark 
N e ff/y/~Nf, to see whether removing the bin would be an improvement. These 
actions mitigate the possibility that the algorithm converges to a false maximum 
with relatively high background. 

• The found set of thresholds are applied to the alternative sets of triggers which 
were kept aside in order to re-estimate unbiased values for N e ff and N b . These 
two numbers, divided respectively by the number of injected events (4320) and of 
time-shifted configurations (400), give the estimates for the efficiency of detection 
and for the average background counts. 



Efficiency and confidence 

• For each shifted configuration of the alternative set, the number of coincident events 
is used to compute the corresponding upper limit at 95% confidence. 

• This entire procedure is repeated for a different couple of detectors and/or a different 
set of MDC injections. 

• After having decided the level of confidence for the eventual rejection of the Null 
hypothesis, we unveil the coincidence counts at zero-delay, and compute for them 
the confidence interval. 

Appendix B. Details about the injected signals 

In this study we assume that the waveforms, in the TT gauge are of the form 




h, 



1 + WP gw T* 



-t/r 



nfgwT\ r(l +e-V(2/^r)) 

cos 2ip — sin 2ip \ 
sin 2i[) cos 2ip J 

W ' B(/ J JcoM^r./;,,,./) 



x 



X 



2 Af gw > 

cos i 0(t) sin (2nf gw t) 



(B.l) 



where the angle ip is an arbitrary polarization and the angle i an arbitrary 
inclination of the angular momentum of the system which originates the burst with 
respect to the line of sight. 
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Notice the different for the h + and h x terms; the reason for multiplying 
cos (2nf gw t) by — |j— ) is to avoid a discontinuity at the beginning of the waveform, 
which would result into an infinite energy, even though h rss would remain finite. 

The polarization ip is uniformly distributed in [0, 2ir), while cosi is uniformly 
distributed in [—1, 1); these choices correspond to assume a random orientation in space 
of the axis of symmetry of the emitting system. 



Appendix B.l. Signal normalization 

The signal normalization is done requiring equation (B.2) for — 0, i — 0: 
h 2 rss = | o °° (\h + (t)\ 2 + \h x (t)\ 2 ) dt . 



(B.2) 



It can also be useful to relate h rss and the energy emitted E assuming a source 
located at a distance r. To this end, we recall the standard definition of the energy flux 



7TC 



f 2 ( ~h + (f) + h,(f) 



dE 

dAdf 2G N 
where dA = r 2 dVL. It is straightforward to compute 



(B.3) 
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B.4) 



16/% 4 r 4 + (4/2 w 7T 2 r 2 + l) 2 + P (8tt 2 t 2 - 32/ 2 w vr 4 r 4 )) 
and then, after performing the integral over frequencies and the angles, one obtains: 
^r 2 h 2 rss (4/>V 4 1 
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2G 
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7 + he 21 f" 
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i.e.: 
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(B.6) 



_io- 2i /v^. 

which means that observing an h rss ~ 10 -21 / \/ H ' z with signals at about 1kHz 
corresponds to a source emitting a fraction ~ 10~ 6 of a solar mass in gravitational 
waves at a distance of lOkpc. 
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